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We have performed a quantitative analysis of recent resonant x-ray scattering experiments carried 
out in the antiferromagnetic phase of V2O3 by means of numerical ab-initio simulations. In order to 
treat magnetic effects, we have developed a method based on multiple scattering theory (MST) and a 
relativistic extension of the Schrodinger Equation, thereby working with the usual non relativistic set 
of quantum numbers l,m,a for angular and spin momenta. Electric dipole-dipole (El-El), dipole- 
quadrupole (E1-E2) and quadrupole-quadrupole (E2-E2) transition were considered altogether. We 
obtain satisfactory agreement with experiments, both in energy and azimuthal scans. All the main 
features of the V K edge Bragg- forbidden reflections with h + k + l =odd can be interpreted in terms 
of the antiferromagnetic ordering only, ie, they are of magnetic origin. In particular the ab-initio 
simulation of the energy scan around the (l,l,l)-monoclinic reflection excludes the possibility of any 
symmetry reduction due to a time-reversal breaking induced by orbital ordering. 



PACS numbers: 78.70.Ck, 71.30. +h 



I. INTRODUCTION 



In the last fifteen years, after the discovery of high- 
temperature superconductivity in cuprates, there has 
been an upsurge of renewed interest in the electronic 
properties of transition metal oxides. 1 Among the differ- 
ent techniques applied to investigate such materials, res- 
onant x-ray scattering (RXS) has proved to be a power- 
ful tool to extract direct informations about magnetic 2 ' 3 
and electronic 4,5 distributions: no other methods have 
the same flexibility to detect lattice, orbital and mag- 
netic anisotropies within the same experimental setup. 6 
Yet, a general theoretical comprehension of some impor- 
tant implications of this technique is still lacking, and 
this circumstance has sometimes led to incorrect conclu- 
sions about the origin of the anomalous signals. 7-9 This is 
particularly true for magnetic RXS, where the absence of 
numerical ab-initio simulations has strongly limited the 
possibility of quantitative investigations. Only recently 
few papers have appeared dealing with this matter. In 
particular Takahashi et al. have used the ab initio lo- 
cal density approximation (LDA) + U scheme, taking 
also into account spin-orbit interaction, to describe such 
phenomena. In this way they were able to calculate the 
magnetic RXS in the energy region of the Ap conduction 
band in KCuF 3 . 

One of the purposes of this paper is to present an alter- 
native method to deal with magnetic phenomena, which 
we believe physically more intuitive, based on a relativis- 
tic extension of the Schrodinger Equation and the multi- 
ple scattering theory (MST). This extension is obtained 
by eliminating the small component of the relativistic 
wave-function in the Dirac Equation and working only 
with the upper component, 12 thereby using the more in- 



tuitive non relativistic set of quantum numbers (l,m,cr) 
for the angular and spin momenta. In order to take into 
account spin-polarized potentials and spin-orbit interac- 
tion in the framework of MST we solve a two-channel 
problem, in each atomic sphere, for the two spin com- 
ponents of the wave-function coupled by the spin-orbit 
interaction. The condition m + a = ml + a' must be sat- 
isfied, due to the local conservation of the z-projection 
of total angular momentum j = I + a. We have imple- 
mented the computer code to calculate magnetic RXS 
into the FDMNES package. 13,14 

Secondly and more important, we present here the first 
application of this method to deal with the case of V2O3, 
which has attracted considerable attention in the past 
four years, 8,9 ' 15-21 due to the peculiar interplay of mag- 
netic and orbital ordering. By performing a quantitative 
analysis of the recent RXS experiments carried out by 
Paolasini et al. 16,17 at the vanadium K edge, we shall try 
to establish the magnetic space group of the monoclinic 
phase, in the endeavor to resolve some controversies that 
are still going on in the scientific community. 

A brief history of the most recent findings about V2O3 
will motivate our work. This compound is a Mott- 
Hubbard system showing a metal-insulator transition at 
around 150 K from a paramagnetic metallic (PM) to an 
anti-ferromagnetic insulating (AFI) phase due to the in- 
terplay between band formation and electron Coulomb 
correlation. 22 A structural phase transition takes place, 
at the same temperature, from corundum to monoclinic 
crystal class. In recent years, on the basis of an old the- 
oretical model by Castellani et al., 23 Fabrizio et al. 15 
suggested that a direct observation of orbital ordering 
(00) could have been possible by means of RXS. Soon 
after Paolasini et al. 16 interpreted the forbidden (111) 
monoclinic Bragg-forbidden reflection as an evidence of 
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such 00. All this gave rise to an intense debate 8 ' 9,17-21 
whose main achievement was to prove the incorrectness 
of the old Castellani et al. model. In fact, non-resonant 
x-ray magnetic scattering experiments 16 and absorption 
linear dichroism 24 showed that the atomic spin on vana- 
dium ions is S=l, and not S=l/2, as supposed in Ref. 
[ 23]. Yet, none of the attempts to explain the origin 
of the (111) monoclinic reflection 8 ' 9,17-20 can be consid- 
ered satisfactory, for reasons that will be clearer in the 
following. Here we just recall that many different physi- 
cal mechanisms have been proposed, from the antiferro- 
quadrupolar ordering of the 3d orbitals, 8 to the combined 
action of the ordered orbital and magnetic degrees of 
freedom; 9 from an orbital ordering associated with a re- 
duction of the magnetic symmetry, 18 to the anisotropies 
in the magnetic octupolar 19 or toroidal distribution. 20 

In this paper we shall focus on a series of 
Bragg-forbiddcn reflections measured by Paolasini and 
collaborators. 16,17 Of particular interest, from a theoret- 
ical point of view, is the set of data of Ref. [17], where 
the different reflections have been collected with partic- 
ular care on their relative intensities, and energy and az- 
imuthal scans have been performed in both the rotated 
(air) and unrotated (crer) channels. Thanks to the wide 
information content of these data, we are able to demon- 
strate that the signal cannot be due to any kind of orbital 
ordering or charge anisotropy and it must be magnetic. 
Moreover we can show that both mechanisms suggested 
in Refs. [ 19, 20] can be at work for the (1,1,1) reflection, 
as there is no extinction rule for either dipole-quadrupole 
or quadrupole-quadrupole transitions. Which of the two 
contributions dominates depends strongly on the kind of 
reflection and on the azimuthal angle. To reach this con- 
clusion ab initio calculations have proved necessary. 

The next section is devoted to the description of our 
method to deal with magnetic phenomena without resort- 
ing to the Dirac equation. We shall work in the frame- 
work of MST with muffin-tin approximation or within 
the finite difference method (FDM), 13 ie, without ap- 
proximation on the geometrical shape of the potential. 
In Sec. Ill we introduce the crystal and electronic prop- 
erties of V2O3 and compare our results with the exper- 
imental data, discussing the physical mechanism behind 
each reflection. Finally we draw some conclusions on the 
possibility of gaining information about orbital and mag- 
netic degrees of freedom by means of RXS. 



II. METHOD OF CALCULATION 



A. RXS and related spectroscopies 



M ng = {lp n \d\ip g } 



(1) 



Here ip g is the ground state and ip n the photo-excited 
state. In the x-ray regime, the operator O is usually writ- 
ten through the multipolar expansion of the photon field 
up to the electric quadrupole term 25 : 



di(o) = gt(o) . ? ^ _ h%i{o) . 



(2) 



where r is the electron position measured from the ab- 
sorbing ion, e^ ) is the polarization of the incoming (out- 
going) photon and fe'(°) its corresponding wave vector. 
In the following dipole and quadrupole will refer to elec- 
tric dipole and electric quadrupole. In RXS the global 
process of photon absorption, virtual excitation of the 
photoelectron, and subsequent decay with re-emission of 
a photon, is coherent throughout the crystal. Such a 
coherence gives rise to the usual Bragg diffraction condi- 
tion, that can be expressed, at resonance, as: 



F= j2e^(f 0a +f' a +ia 



(3) 



where R a stands for the position of the scattering ion a, Q 
is the diffraction vector and /o is the usual Thomson fac- 
tor. The resonant part, f' + if", is the anomalous atomic 
scattering factor (ASF), given by the expression 26 : 
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(4) 



(E n - Eg) - i l ~f 

Here Two is the photon energy, m e the electron mass, E g 
the ground state energy, and E n and T n are the energy 
and inverse lifetime of the excited states. In practice, the 
intermediate states tp n are in the continuum, normalized 
to one state per Rydberg. For this reason it is useful to 
label them by their energy E and re-express Eq. (4) as: 



/' + if" = -m e uj 



2 E 



M°*M % 

ng ng 



e f E — E„ — Tiuo 



dE (5) 



using the fact that, at resonance, E n — E g ~ Two. Ep is 
the Fermi energy. In Eq. (5) the summation over n is 
now limited to states having the energy E. 

It is now straightforward to make the connection with 
another spectroscopic technique: x-ray near edge ab- 
sorption (XANES). In this case the cross-section sim- 
ply corresponds to the imaginary part of the ASF when 
o — Oi (forward scattering). Choosing prefactors in 
order to have /" in units of the classical electron radius, 
r ~ 2.82 • 10~ 15 m, and the absorption cross section a 
in megabarn, we have the relation: 



Core resonant spectroscopies are described by the vir- 
tual processes that promote a core electron to some 
empty energy levels. They all depend on the transition 
matrix elements of the operator O expressing the inter- 
action of electromagnetic radiation with matter: 



-4tt 10 



22 a a"m e c- 



fiuj 



f" 



(6) 



where ao is the Bohr radius, a is the fine structure con- 
stant and c is the speed of light. The close connection 
between the two spectroscopies is thus evident. 
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The use of RXS can have some advantages in exploring 
electronic properties compared to the absorption tech- 
niques. First, with RXS it is possible to select different 
relative conditions in incoming and outgoing polariza- 
tions and wave vectors, and this gives more opportunities 
to probe the magnetic and electronic anisotropics of the 
material. Secondly, RXS can be more site selective than 
XANES, due to the Bragg factor (3). For instance, it 
is possible, in magnetite, to be sensitive only to octahe- 
dral Fe 3+ -sites and not to tctrahedral Fe 35+ -sites 27 or, 
in manganites, to probe only Mn 3+ -ions and not Mn 4+ - 
ions. 28 This can be achieved at reflections that are Bragg- 
forbiddcn off-resonance, but become detectable, in reso- 
nant conditions, because of non-symmorphic symmetry 
elements (glide planes, screw axes) in the crystal space 
group. 29 In these cases the Thomson factors /o (Eq. 3) 
drop out and one is just sensitive to the changes of the 
ASF (Eq. (4)) that depend on the relative electronic 
and magnetic anisotropies on the ions related by such 
non-symmorphic symmetry elements. It is indeed in such 
kinds of reflections that orbital ordering or magnetic scat- 
tering were reported and they will be the main subject 
of the present study. 

B. Cartesian Tensor approach 

In Sec. Ill we shall deal with the symmetry operations 
of the crystal space group of V2O3, in order to evaluate 
the structure factor for RXS. Because of this, it is very 
useful to write explicitly the dipole and quadrupole com- 
ponents of the matrix elements in Eq. (5). Remembering 
that incoming and outgoing x-rays can have different po- 
larizations and wave vectors, we get: 

n afi 

a/37 

+ \ E (7) 

where a, (3, 7 and 5 are cartesian coordinate labels and 
D a p, Ia/3-y and Qap-ySi the dipole-dipole (dd), dipolc- 
quadrupole (dq) and quadrupole-quadrupole (qq) contri- 
butions, respectively. Their explicit expression is given 
by: 

D a p = ^{lpg\rMn){tpn\r0\^g) 
n 

E (<Pg\r a \ipn)(ipn\rpr 7 \ip g ) 

n 

= ^2('<l>g\r a rp\Tpn}(ipn\r~ ( r s \iP g } (8) 

n 

Note that at a K-edge, in absence of spin-orbit interac- 
tion and time- reversal breaking potentials, D, I and Q 



are all real. 



C. Calculation of the excited states 

The main difficulty in the ab-initio evaluation of the 
ASF is the determination of the excited states. Two dif- 
ferent procedures, MST and FDM, have already been 
developed for non-magnetic cases, and the corresponding 
packages 14 used in some XANES 30 and RXS 10 ' 31 appli- 
cations. In the following we introduce a relativistic ex- 
tension of the ab initio cluster approach, for both MST 
and FDM, that is still mono-electronic but includes the 
spin-orbit interaction and allows to handle magnetic pro- 
cesses. The importance of the FDM procedure to avoid 
the muffin-tin approximation has already been reported 13 
and will not be detailed here again. Note only that with 
FDM the intermediate states are calculated by solving 
the Schrodingcr equation without any approximation on 
the geometrical shape of the potential, ie, we do not need 
to use, as usually in MST approaches, a spherically av- 
eraged potential in the atomic spheres and a constant 
among them. This point is essential to study orbitally 
ordered materials, where an anisotropic orbital distribu- 
tion is the key-feature of the system and no spherically 
averaged potentials would be adequate for its description. 
Unfortunately, large cluster calculations with FDM are 
prohibitively long in computing time and require a lot of 
computer memory: that's why, when appropriate, MST 
is preferred. 

In order to treat magnetic effects in MST we have con- 
sidered the relativistic extension of Schrodingcr equation 
obtained by solving Dirac equation exactly for the upper 
component of the wave-function, as described by Wood 
and Boring. 12 This approach ensures that the singular- 
ity of the potential at the nuclear center is of centrifugal 
type (<~ r~ 2 ) even in the presence of spin-orbit terms 
(so that the usual fuchsian method of solution can be 
used) . It also allows us to work with the usual non rela- 
tivistic set of quantum numbers (I, m, a) of angular and 
spin momenta, which we find much more intuitive than 
the relativistic set. Our one particle basis will therefore 
comprise spin-orbit coupled core states of the form: 

*k;^ = WLm*;Jc) (9) 

= R c 3 (r) Y lc m c X« (l c m c l/2a\j,j z ) 

m c <7 

and valence excited states which are multiple scattering 
solutions of the Schrodinger equation with spin-polarized 
potentials. In atomic units: 

{v 2 + k 2 e - V (r) - Vitfs, - W 2 {r)t s) * fe - s (f) = 

(10) 

These latters are supplemented by the scattered wave 
boundary conditions: 
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*E e ,a(0 = ~ fike.s-r.s'f—xs' (11) 

In Eq. (9) (l c m c l/2a\j,j z ) are the Clebsch-Gordan co- 
efficients. In Eq. (10) Vo(f) is the average of spin up and 
spin down potentials, V\(r) their difference and V2(r) the 
usual spin-orbit term. Such potentials already embody 
the relativistic corrections due to the reduction of the up- 
per Dirac component of the wave function to an effective 
Schrodinger equation. 12 

In the muffin-tin approximation the solution inside the 
i-th atomic muffin-tin sphere can be written as: 

*£.,,(*) = EE E Rt/(n)B$ m , a ,(k e , s )Y lm (n) Xa 

l ma m'a' (12) 

where at the muffin-tin radius pi the radial functions 
R\ma ( r i) match smoothly to the following combination 
of Bessel and Hankel functions via the atomic ti m<7 j m > a > 
matrices (defined below): 

R\t:'^) = jl{kpiKL,lm>a> ~ iSlma,lmia>ht{k Pi ) (13) 

With a proper normalization of these radial functions 
to one state per Rydberg, the scattering amplitudes 
B LrT (k e , s) obey the equations (writing L for lm): 

E E M t,L>a'B j L ,Ake, S) = 8 s jY L (kY k ^ 
j L'a' 

where 

M La,L'a' = {^ImaM m< a' &i]&W &m+ajn' +a> + G % [ \ l ,5<j,<ti} 

is the usual multiple scattering matrix, generalized to 
spin variables, and Ri denotes the position of the i-th 
atom in the cluster with respect to the origin of the co- 
ordinates. The atomic i-matrix U ma j m i a i describes the 
scattering amplitude of an electron impinging the atomic 
potential with angular momentum I, azimuthal compo- 
nent m and spin a into a state with quantum numbers 
l,m',a'. The conservation of the total angular momen- 
tum j ' = I + s and its z-projection implies the constraint 
m + a = m' + a'. Note that also I is unchanged in the 
scattering process, since in the muffin-tin approximation 
the potential has spherical symmetry. 

By introducing as usual 33 the scattering path operator 
T La L'a' as the inverse of M % [ a U(jl , the solution for the 
scattering amplitudes B La is given by: 

Bi„{k e ,8) = 5 as Y J Y. T t,L'/Y L '{ky k ^ (14) 

3 L' 

The scattering path operator Tj? L , a i represents the 
probability amplitude for the excited photoelectron to 
propagate from site i, with angular momentum L and 
spin a, to site j with angular momentum L' and spin 
a'. It is the obvious generalization of the corresponding 
spin-independent quantity. 33 



We now have all the ingredients to perform the inter- 
mediate sum in Eq. (5). Using the expression given Eq. 
(12) for the state \ip n ) in the matrix elements of the nu- 
merator, the sum over n becomes an integral over the 
escape direction of the photoelectron, fc e , at fixed energy, 
E = fc 2 . When we perform the integral over the energy, 
we get the product of two smoothly varying radial matrix 
elements of the type (R\™£ {r)\r h '\Rj(r)) times the ap- 
propriate angular (Gaunt) coefficient. The exponent v is 
determined by the polar order of the transition. This is 
further multiplied by the angular integral of the product 
of the two scattering amplitudes B l La , B J L , a , centered on 
the absorbing sites. This latter can be simplified by the 
use of a generalized optical theorem 32 : 

]T f dk e Bi a (k e ,s)Bi >a ,(k e ,s) = 3r^ iv , (15) 

As a consequence, the knowledge of r at all relevant en- 
ergies, and of the radial solutions R^J (r) on the pho- 
toabsorber are sufficient to calculate the ASF. 

In the case of V2O3, in order to construct the spin- 
polarized atomic potentials, we have used the prescrip- 
tion by von Barth and Hedin : they were derived from 
the non-self-consistcnt spin-polarized charge density ob- 
tained by superimposing the atomic charge densities with 
two magnetic electrons on each Vanadium ion, as sug- 
gested by the experimental data. 16 ' 24 We believe that this 
approximation is not far from reality, as shown a posteri- 
ori by the goodness of the results. Given the potential, it 
is straightforward to calculate the atomic spin dependent 
t-matrix: we just solve the two-channel problem arising 
from the Schrodinger equation (10) due to the presence of 
the spin-orbit potential. In this way this latter is treated 
on the same footing as the other potentials and not in 
perturbation theory as usually done. 

The relativistic extension in the FDM scheme follows 
closely the previous treatment for MST. The space is par- 
tioned, as in the multiple scattering approach, in three 
regions 14 : i) An outer sphere surrounding the cluster 
where one impose the scattering behaviour of the wave 
function given in Eq. (11). ii) An atomic region made 
up of spheres, around each atom, with radii much smaller 
than the muffin-tin radii (of the order of 1 a.u. or less). 
Here the charge density is to a very good approximation 
spherically symmetric, due to the presence of the core 
electrons. Hi) Finally, an interstitial region where the 
Laplacian operator of the Schrodinger equation (10) is 
discretized and the solution is generated on a greed with- 
out any approximation on the geometrical shape of the 
potential. By imposing a smooth continuity of the over- 
all wave function across the boundaries of the three re- 
gions, one can determine both the expansion amplitudes 
of the wave function inside each atom and the scatter- 
ing T-matrix of the whole cluster in Eq. (II). 13 Note 
that in this case the spin-orbit potential is not spheri- 
cally symmetric. We have checked that for a small cluster 
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with close-packed geometry, this approach and the mul- 
tiple scattering approach with muffin-tin approximation 
provide almost identical cross sections. Since V2O3 is a 
close-packed structure, if one neglects vanadium voids, 
we have used MST in the muffin-tin approximation for 
the magnetic calculations. The FDM scheme was instead 
used in the calculation of the ASF in the case of orbital 
ordering. 



III. THE CASE OF V 2 3 

In this section we specialize to the case of V2O3, to in- 
vestigate whether the RXS reflections observed by Pao- 
lasini and collaborators 16 ' 17 are of magnetic or orbital 
nature. Our result is that magnetic ordering can explain 
the experimental data with a reasonably good agreement 
(even if not perfect, for reasons that will be discussed 
below). On the other hand, we are able to demonstrate 
that any 00 origin of the signal has to be excluded, as 
it would not fit the energy scan at dipolar energies. 




FIG. 1. Monoclinic cell in the AFI phase of V2O3. The 
numbers label the eight vanadium ions of the unit cell. The 
arrows show the spin orientations. a m , b m and c m are the 
monoclinic axis. Only V ions are represented: white and gray 
circles indicate that they are surrounded by differently ori- 
ented oxygen octahedra. Two basal planes of the hexagonal 
corundum cell are also shown. 

Next subsection is devoted to the description of the 
magnetic space group in the AFI phase, as well as the 
derivation of the anomalous RXS amplitude. In subsec- 
tion B we analyze the ab-initio calculations for the mag- 
netic RXS, through the procedure previously developed. 
Finally, in the last part we show that the consequences 
of a time-reversal breaking 00 are not compatible with 
the experimental data, thus putting a severe constraint 
on the possible 00 theories for V2O3. 



A. Crystal and magnetic structure 

As a first step we shall try to establish the magnetic 
symmetry group in the AFI phase. Our starting point 
is the crystallographic structure, as given by Dernier 
and Marezio 36 together with the magnetic data given by 
Moon 38 and later confirmed by Wei Bao et al. 39 

Referring to the frame and the numbering of vanadium 
atoms of Fig. 1, we can divide the eight atoms of the 
monoclinic unit cell into two groups of four, V\, V2, V3, V4 
and Vi', V2', V31 , V4', with opposite orientations of the 
magnetic moment. The two groups are related by a body- 
centered translation. Vanadium-ion magnetic moments, 
indicated by arrows in the figure, lie perpendicular to the 
6 m -axis at an angle of 138° away from the monoclinic c m - 
axis. Neglecting the magnetic moments, the previous two 
groups of four atoms with their oxygen environments are 
translationally equivalent. 

We can infer from these data that the magnetic space 
group can be written as P2/a + T{E\t }P2/a, where T 
is the time-reversal operator and to the body-centered 
translation. The monoclinic group P2/a contains four 
symmetry operations: the identity E, the inversion j, 
the two-fold rotation about the monoclinic b m axis C2b 
and the reflection mj, with respect to the plane perpen- 
dicular to this axis. All these operations are associated 
to the appropriate translations, as shown in the following 
tabic: 



atom 


position 


spin 


sym. 


translation 


Vi 


(u,v,w) 


T 


E 





v 2 


(-u,-v,-w) 


T 


i 





V3 


(| + u, —v, w) 


T 




(|,o,o) 


v 4 


(5 -u,v,-w) 


T 


fc 2b 


(|,o,o) 


Vv 

Vy 


(i + U, ~ + V, i + W) 

(1 - u > \ - v >\ - w ) 

(u, 5 -v, ! + «;) 


1 
I 
I 


f 

fi 

m b 


(1 1 1) 

V 2 ' 2 ' 2 1 

V 2 ' 2 ' 2 / 

(o,|,|) 




(-U, \ + v, \ - w) 


I 


C2b 


(o,|,|) 



TABLE I. 

Here u = 0.3438, v = 0.0008, w = 0.2991 are the frac- 
tional coordinates of the atoms in unit of the monoclinic 
axis. Notice that it is sufficient to calculate the ASF for 
atom number 1, since all the others can be deduced by 
applying the symmetry operation indicated in the fourth 
column of the table. 

In the following, we re-analyze the extinction rules for 
the different reflections (h, k, I) measured by Paolasini et 
al. 16 ' 17 in terms of the crystal tensors in the dd, dq and 
qq channels (Eq. (8)). All the recorded reflections have 
the same incoming polarization, a (electric field perpen- 
dicular to the diffraction plane), whereas the outgoing 
polarization is analyzed both in the a and it channels 
(this latter with the electric field in the diffraction plane). 
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The azimuthal scans are then registered by rotating the 
sample around the diffraction vector. We refer to Pao- 
lasini et al. [ 17] for the definition of the azimuthal origin: 
the situation where the scattering plane contains a and c 
hexagonal axes corresponds to an azimuthal angle of 90 
degrees. 

From equation (3) and Table I, neglecting the very 
small component of R along b m , (v = 0.0008), we get for 
each reflection (h, k, I) the expression: 



hkl 



(l + (-)''Tm fc )(e ,fe! +- 

(i + H 



h+k+l 



</) 



(16) 



where R ~ (u, 0, w) is the position of the V\ atom, 
4>hkl = Q R = 2ir(hu + Iw) and /i stands for the ASF as 
defined in (Eq. 5). Note that /i is a scalar and we adopt 
the convention, when we say that a symmetry operator 
acts on /i, that it acts only on the tensor components 
(dd, dq and qq) given by Eq. (8). Focusing on the re- 
flections with h + k + l =odd and noticing that D a p and 
Qa/3-yS are inversion-even, while I a p 7 is inversion-odd, we 
find that only the following tensor components contribute 
to the signal: 

= 4*(1 + (-i) h+n y +1 )cos^ hkl )Q(D aP ) 
:$ 7 = 4(1 + sin(^ fei )3(/^ 7 ) 



Q^U = 4i(l + COs(0 fcfcI )9(Oa/J 7 «) ( 



17) 



where is the number of j/ labels among the tensor in- 
dices a, /?, 7, and <5. Thus, in order to have a signal, 
h and n y must have different parity. Notice that all the 
quantities in Eq. (17) arc magnetic, as only imaginary 
parts of cartesian tensors are involved. Moreover, the 
three scattering amplitudes are all purely imaginary, as 
the dq polarizations carry an extra imaginary unit (see 
Eq. 7). As a consequence, they all interfere. A sepa- 
rate analysis of the three tensors of Eq. (17) gives the 
following indications. 

In the case of dd tensors, when h is odd, the only 

When h is 



•\hkl 



non-zero component is 2?£* — T> 
even the non-zero components are T> 
and X>$« - Z£« ex (L x ). 



(Ly) 

hkl _ 

xy 



V™ cx (L z ) 
As the magnetic moment 



direction is perpendicular to the frm-axis, \Ly) is zero. 
Thus for odd h no signal is expected in the dipolar re- 
gion. On the contrary, when h is even, a dd contribu- 
tion is present. Of course, no a a magnetic scattering 
is allowed. 26 These facts explain why the experimental 
spectrum for the (2, 2, shows structures at the 4p 
edge, contrary to the (2, 2, X) aa reflection and to the 
(1, 1, 1) and (3, 1, 1) reflections for both polarization con- 
ditions. Notice that these results had been already de- 
rived with similar methods. 18-20 

These remarks do not apply for qq and dq tensors. 
They always contribute, for all the investigated reflec- 
tions, in both a a and air channels, except at some very 
specific azimuthal angles. At the K edge these tensors 



measure the transitions to the states with pure 3d or hy- 
bridized 3d-4p character. For this reason they have a 
finite value just close to the Fermi energy. The physical 
quantities measured by these operators can be identified 
as the octupolar magnetic moment for the qq term and 
the toroidal or quadrupolar magnetic moment for the dq 
tensor. 

B. Analysis of the magnetic signal 

In order to get the absolute intensities and shape of 
the spectra, we need now to resort to ab initio calcu- 
lations. We performed such calculations for the crystal 
and magnetic structures given in Table I. Thus, we do 
not neglect the small v value but we shall see that the 
conclusion given in the previous subsection is not modi- 
fied. The potential is calculated using a superposition of 
atomic densities obtained from an atomic, self-consistent 
Hartree-Fock calculation, with a 3d 2 , spin 1, configura- 
tion. The MST approach is used with different cluster 
radii from 3.0 up to 7.2 A, ie, from a VOq molecule to a 
cluster containing 153 atoms. 
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FIG. 2. Comparison of the experimental (dots) and cal- 
culated (continuous line) spectra for the XANES and the 
(2, 2, l) CTVr reflection at the vanadium K edge. Calculations 
are performed with three different cluster radii (3, 4.3 and 7.2 
A). The pre-edge feature is also shown before the convolution 
for the 4.3 A case. The first structure, indicated by an ar- 
row, is at the energy where the 3d feature of the Bragg peaks 
appears. 

We find that, in order to get the shoulder at 5475 eV 
in both XANES and (2, 2, l) a7T spectra, as well as the 
5487 eV shoulder in the (2, 2, 1)^, we need the biggest 
cluster radius (7.2 A - see Fig. 2). Nevertheless the main 
features are present for all Bragg peaks even for the VOq 
molecule calculation. 
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The behaviour of the azimuthal scans against the clus- 
ter radius looks more complex. The agreement improves 
up to the 4.3 A, corresponding to 33 atoms, and then de- 
creases when the cluster radius is further increased. The 
reason for this is discussed below. For the moment we 
keep the 4.3 A radius and compare such energy and az- 
imuthal spectra with the experimental ones, as shown in 
Figs. 3 and 4. 
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FIG. 3. Comparison of the experimental (dots) and calcu- 
lated (continuous line) energy scans for three different reflec- 
tions. On top is also shown the XANES spectrum. The calcu- 
lation is performed on a 4.3 A cluster radius and includes the 
magnetic ordering. The intensity unit is the electron number 
squared. Note that the experimental spectra were recorded 17 
with a special care on the relative amplitudes. For the (2, 2, 1) 
and (3, 1, 1) reflections the azimuthal angle as denned by Pao- 
lasini and coworkers is ip — 15°. All the mean features are 
present. An extra scale factor is applied on the experimental 
dots of the (1, 1, l) at7 reflection. 



Looking at the spectra shown in Fig. 3, we can claim 
a satisfactory experiment-theory agreement for the dif- 
ferent reflections. In particular we get the very thin res- 
onance in the 3d — i 2 g energy range at the (1,1,1) and 



(3, 1, 1) peaks for both aa and air polarization conditions 
and for the (2, 2, 1) (J(T reflection. The broad intensity in 
the 4p energy range in the (2, 2, l) CT7r is also reproduced. 




0.5 - 



& 1 
a 
o l 



1 - 



1 1 ' 1 1 1 1 1 1 1 1 1 1 1 1 ' I ' 1 ' 1 1 1 ' 1 ' I ' 1 ' 1 1 1 

(1,1,1)071 



I | I I I I | I I I T | T I I I | I I I tp 1 1 I | I I I TffT 
x0.5 (l,l,l)oo 




(2,-2, l)oo 




-150-100 -50 50 100 150 

Angle 

FIG. 4. Azimuthal scans of the same three reflections 
shown in Fig. 3. The main features of the experiment (dots) 
are obtained in the calculation (continuous line). The scans 
are performed respectively at 5464 eV for the (1,1,1) and 
(3, 1, 1) reflections and 5465 eV for the (2, 2, 1). The intensity 
unit is the number of electrons squared. The azimuth 90° cor- 
responds to the situation where the scattering plane contains 
the c hexagonal axis. 

The intensity ratio between aa and air channels and 
between the various reflections are also quite good con- 
sidering the difficulties on both experimental and theo- 
retical sides, even if for the (1, 1, reflection an extra 
factor improves the agreement with the experiment. Note 
that, differently from all other reflections, the (2, 2, l) a7r 
can contain some non-resonant component not taken into 
account in our calculation. Moreover, a non-zero offset 
seems to be present at (2, 2, l) aa on the experimental 
side, and this is probably responsible for the discrepancy 
in the corresponding energy and azimuthal scans. The 
azimuthal scans, shown in Fig. 4, are reasonable for all 
air reflections. Only the (1, 1, l) aa is much less satisfy- 
ing: yet, the corresponding data belong to the first ex- 
periment (Ref. 16, while the others belong to Ref. 17), 
and this does not allow to be sure about the relative in- 
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tensity as in the other cases. Finally, we are also able to 
reproduce the it and 2ir periodicity, expected from the 
low symmetry of the compound, for the oo and an chan- 
nels, respectively. 
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FIG. 5. Magnetic intensity of the various reflections at 
the vanadium K edge. The total (continuous line), dd 
(dash-dotted line), dq (dash line) and qq (dotted line) sig- 
nals are shown, at the energy E=5465 eV 

We judge the agreement between theory and experi- 
ment quite satisfactory. Indeed the interference between 
the dd, dq and qq parts makes the problem really tricky. 
Foe example, p orbitals are far more delocalized than d 
ones. For this reason they react in a different way to the 
core hole screening. In a first approximation d states are 
more shifted towards lower energy than p states. Thus, 
the dq terms, which probe the hybridization between 
p and d orbitals, are strongly influenced by the details 
of the electronic structure around the photo-absorbing 
ion. If we consider that V2O3 is a strongly correlated 
electron system, with a quantum-entangled ground-state 
wave function, 9 ' 18 ' 20 we find quite surprising that we have 
reached such a good agreement in the 3d energy region 
with a one-electron calculation. Our idea is that this 
can help understanding which features can be explained 
in terms of an independent particle approach and which 
cannot (see below). Note that the localization of the d 
orbitals also explains why the azimuthal scans look bet- 



ter when the calculations are performed using the smaller 
4.3 A cluster radius. To end this discussion we stress that 
energy spectra profiles around 5464 eV are very sensitive 
to the azimuthal angle and, vice-versa, azimuthal pro- 
files are very sensitive to the energy. This is due to the 
fact that qq and dq components have both strong (and 
different) angular and energy dependence. There are 
some angles where the calculated profiles appear doubled, 
with about 1 eV between the two peaks. For instance, 
this is what happens at the experimental azimuth of the 
(3, 1, l) aa energy spectrum: a simple shoulder appears at 
5465 eV, well described by our calculation. 

Finally, we want to comment about the two papers 
(by Lovesey et al. 19 and by Tanaka. 20 ) that have already 
pointed out the magnetic origin of the h + k + I =odd 
reflections. These works disagree on the physical mech- 
anism of the signal, as for Tanaka only dq terms con- 
tribute, while Lovesey and coworkers attribute the signal 
to a pure qq reflection. Our results on this point show 
that both dq and qq channels can contribute to the global 
intensity depending on the reflections (whether h =odd 
or h =even) and on the azimuthal angle (see Fig. 5). 
At the prepeak energy, E— 5465 eV, the dq term rep- 
resents the strongest contribution for the (1,1,1) and 
(3, 1, 1) reflections, even if the qq term is not completely 
ncgligcablc. For the (2, 2, l) CTO -, on the contrary, the qq 
contribution is clearly the dominating one. Finally, it 
is interesting to note that for the (2, 2, l) av reflection 
the dd term plays the major role even at the 3d en- 
ergies. Anyway the qq contribution has always a non- 
negligeable influence on the total signal. Comparing our 
findings with those of Refs. 19, 20 we can state that there 
are no strict extinction rules regarding dq and qq terms 
for all these reflections. Nonetheless, the main contribu- 
tion to h + k + I =odd, h —odd reflections comes from 
the dq channel (as in Ref. 20), while the main term in 
h + k + I =odd, h =evcn reflections is of qq origin (as in 
Rcf. 19). Such a tensor analysis allows also to identify a 
feature that is not correclty described by our one-electron 
calculations, ie, the direction of the magnetic moment. 
Indeed, a small L y component (about 10% the global mo- 
ment) is found to contribute to the h =odd signal (see 
Fig. 5), contrary to what expected by our previous theo- 
retical analysis. This result can be explained by noticing 
that in order to obtain the correct direction of the mag- 
netic moment a molecular, correlated ground-state wave 
function is needed, 20 which is far beyond the possibilities 
of our monoelectronic approach. 



C. Analysis of the "orbital" signal 

In this subsection, we analyze the effect of a time- 
reversal breaking OO on the signal, in order to deter- 
mine whether it can affect our previous results. This 
is not a secondary issue since the original experimental 
interpretation suggested that the (lll)-monoclinic reflec- 
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tion was due to the OO 16 and a big controversy around 
this point arose in the literature. 8,9,18,20 The reason why 
we focus on a "time- reversal breaking" OO is that, for 
h + k + I =odd, the signal is proportional to (E — T)/i, 
as clear from Eq. (16), so that a non- magnetic signal 
is allowed only when the time-reversal symmetry is bro- 
ken. This is obtained, for example, when Vj and W 
have a different orbital occupancy, 18 and we call such a 
situation time-reversal breaking OO. Different types of 
orbitally-ordered ground states have been suggested in 
the literature, 9,18,23 but none of them possesses such a 
feature. As a consequence, non-magnetic signals at the 
(lll) m reflection are not expected, unless some ad hoc 
hypoteses are made, as recognized in the discussion of 
Sec. VII. D of Ref. [ 18]. One possible way to get a signal 
of OO origin is to consider the first excited state found in 
Ref. [ 18], which has a magnetoelectric (ME) symmetry. 
It is very close to the ground state (~ 1 meV) and, be- 
cause of this, it could be partly occupied. Independently 
of this, if we consider the two maximal ME subgroups of 
the full space groups, ie, P2'/a (with symmetry elements 
E, TI, TC26, rhb) and P2/a' (with symmetry elements 
E, TI, (726, Trhb), it is in principle possible to get a non- 
magnetic signal at the (lll) m reflection. In fact, in both 
cases there are two subgroups of four ions whose elec- 
tronic densities have anisotropics not connected by time- 
reversal nor any other symmetry operation. Consider, 
for example, the case of P2/a', that could be responsible 
for non- reciprocal dichroism 21 : the two groups of atoms 
(Vi, V21, V3, V41) and (Vy, V2, Vy, V4) are independent 
and the structure factor, Eq. 16, becomes 



V2 



J= hkl = (1 + (-) ft Tm t )(e'*"' + (-) 



h+k+l. 



l TI) 



q-hkl 



(h + (-) h+k+l fv) (18) 

The analogous of Eq. 17 for reflections with h + k + 
I =odd is now: 

V™ = 4i(l + (-l) fc+n »)sin(^ fcI )K(Da/j) 
+4i(l + (-l)' l +"* +1 ) cos(^ fcJ )9(Da/j) 
H = 4(1 + (-l) h+n »)cos(0 h «)Si(J a/ , 7 ) 
+4(1 + (-l) h+n y +1 )sm(4> hkl )^{I aPl ) 
Q^ s = 4i(l + (-l) h+n y) sin(0 fcjH )R(Qa/J 7 *) 

4i(l + {-\f+ n v+ 1 ) cos(<j> hkl )Z(Q a0J s) (19) 

Thus, the breakdown of time-reversal symmetry, with 
the introduction of the two ASF /1 and fy, has made 
possible also non-magnetic scattering, through the real 
part of the tensor components. Notice that, due to the 
magnetic space group, both magnetic and non-magnetic 
scattering can in principle interfere at particular polar- 
ization conditions and reflections. 

In order to calculate the signal explicitly, we need to 
express the orbital wave functions of Vi and Vi' ions. 
Starting from the orbital occupancy of Refs. 9, 18, we 
get, in the reference frame of Fig. 1: 



1 



1 



V2 



1 



— ,jd xy + r-dy Z + y^^z 2 g d xz ^d x 2_ y 2 

(20) 



From this orbital occupancy, we can calculate the elec- 
tron density. Then, solving Poisson's equation, we get 
the Coulomb potential. The energy-dependent exchange- 
correlation potential is then obtained with conventional 
procedures. 30 The calculation of the ASF is performed 
for the Vi and Vy atoms with the FDM option. We kept 
a rather small cluster radius (3.0 A), but included the 
oxygen octahedra as well as the first shell of four vana- 
dium neighbors. The ASF of the other atoms are then 
deduced by means of the magnetic symmetries of Table 
I for the two separate subgroups (see also Ref. [ 18], sec. 
VII. D, for a derivation) and the total scattering ampli- 
tude is obtained using Eq. (3). 
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FIG. 6. Intensity of the (1, 1, 1) CTCT (full line) and (1, 1, l) CT7r 
(dotted line) reflections at the vanadium K edge calculated 
with the orbital ordering of Eq. (20). The Bragg peaks have 
the same energy scale as the XANES spectra (full line). The 
big structure at the 4p edge is not present in the experiment. 
The calculation is performed with a 3 A cluster radius. 

Our results for the (1,1, l) m reflection are shown in 
Fig. 6. As in the experiment, 16 a structure is obtained at 
the 3d-energy level. Yet, a much bigger signal is present 
at the 4p edge where the experimental intensity shows no 
features. Such a disagreement is the proof that the OO 
cannot be responsible of the signal: The time-reversal 
breaking OO makes the (lll) m reflection allowed not 
only at the 3d-energy level, but also at the 4p-energies. 
Note that this result remains valid for any time-reversal 
breaking OO, not only the one proposed in Eq. (20), 
as should have been qualitatively expected. Indeed, it 
has already been shown 10,40 that the 3d-4p hybridization, 
even if small, must give a contribution to the dipolar K 
edge signal. In fact, previous numerical simulations in 
another compound, LaMnOa, have shown that the effect 
of a symmetry-breaking OO of 3d orbitals produces a sig- 
nal also at the 4p-energy levels. In this latter case the 
contribution coming from the Jahn- Teller distortion over- 
whelms the OO signal. However in the present case, since 
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Vi and Vi' have the same local distortion, such a signal 
should have been detected, as clear from our numerical 
simulation shown in Fig. 6. Thus, the experimental evi- 
dence of absence of any signal at the 4p energies proves 
two facts: first, there is no time-reversal breaking 00; 
second, the (111) TO cannot be due to 00. A direct conse- 
quence of the previous analysis is that no magnctoclcctric 
s ubgroups are compatible with the (lll) m energy scan, 
as they all break the time-reversal symmetry, in keeping 
with the negative experimental evidence for magnetoelec- 
tricity in V 2 3 . 41 This strongly supports the results of 
Rcf. '[21]. 

Notice that a previous study 42 on the same reflection 
in V2O3 ended up with the conclusion that the 00 was 
responsible of the signal. Unfortunately that calculation 
was performed only around the 3d-cnergy range, where 
the agreement in the azimuthal scan was quite good and 
the big feature at the 4p energies was not detected. The 
lesson to be drawn from this fact is that the azimuthal 
scan around the diffraction vector is not always a fun- 
damental feature in determining the origin of the reflec- 
tion. In fact, it usually reflects more the geometry than 
the electronic properties of the material. This is a rather 
general comment, not limited to V2O3: when there are 
only one or two ab initio independent factors (the tensor 
components), the angular properties arc much more de- 
termined by the geometry (ie, the symmetries) than by 
the dynamics (ie, the relative weight of the radial matrix 
elements). In this situation only the combined analysis 
of both energy and azimuthal scans is a reliable tool to 
investigate the electronic origin of RXS phenomena. 



IV. CONCLUSION 

The main results of the present paper can be inferred 
from the conclusions of the last two subsections dis- 
cussing the results appeared in the literature. In fact, 
up to now, it was possible to classify the interpretations 
of the forbidden Bragg reflections with h + k + I =odd 
and h =odd along two main lines of thought: those who 
explained such reflections in terms of orbital ordering, 
both associated 18 or not 9 with a reduction of the mag- 
netic symmetry, and those who were inclined to a mag- 
netic origin. 19,20 A third explanation, that of the anti- 
ferroquadrupolar ordering proposed in Ref. [ 8] had al- 
ready been ruled out in Ref. [ 18]. Whether the for- 
mer or the latter interpretations were correct, was not 
simple to decide on the basis of the azimuthal scans, 
only. In fact, such scans, for a given multipolar chan- 
nel, measure the crystal symmetry, rather than the elec- 
tronic origin of the reflections. The proof of such a 
statement is that three different mechanisms (magnetic 
in the qq channel 19 , magnetic in the dq channel 20 and 
non-magnetic, via orbital ordering 42 ) gave all a rather 
good agreement in the azimuthal scan at the pre-K edge. 
What definitively rules out the 00 origin of such reflec- 



tions is the energy scan shown in Fig. 6. As noted previ- 
ously, a direct consequence of this result is that no reduc- 
tion of the magnetic space group P2/a + T{E\t }P2/a 
is present. 21 ' 41 

The second important result of the present paper lies 
in the fact that we were able to make a complete ab-initio 
analysis of all h + k + I —odd reflections, starting from 
the crystal and magnetic structure, only. We showed that 
their origin (whether dd, dq or qq) strongly depends on 
the kind of reflection (ie, h =even or h =odd), photon 
energy and azimuthal angle. For h =odd, the dq channel 
is predominant, in keeping with the cluster calculations 
of Tanaka. 20 In spite of this particular agreement, we be- 
lieve that it is important to stress that our approach goes 
beyond the cluster calculations performed in Ref. [ 20] 
as well as the simple fitting procedure of Ref. [ 19]. In 
fact only by means of an ab-initio procedure it is possi- 
ble to cover all the main experimental evidence, ie, the 
azimuthal behaviour, the energy profiles, and the order 
of magnitude of the intensity, that is well in keeping with 
the rough estimate given in Ref. [16]. 

We would like to acknowledge L. Paolasini for several 
useful discussions. 
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